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Abstract 



For monochromatic solutions of D'Alembert's wave equation and 
Maxwell's equations, we obtain sharp bounds on the sup norm as a 
function of the far field energy. The extremizer in the scalar case is 
radial. In the case of Maxwell's equation, the electric field maximizing 
the value at the origin follows longitude lines on the sphere at infinity. 
In dimension d = 3 the highest electric field for Maxwell's equation is 
smaller by a factor 2/3 than the highest corresponding scalar waves. 

The highest electric field densities on the balls -Br(O) occur as 
R — > 0. The density dips to half max at R approximately equal to 
one third the wavelength. The extremizing fields are identical to those 
that attain the maximum field intensity at the origin. 

Key words: Maxwell equations, focusing, energy density, extreme light 
initiative. 
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1 Introduction. 



The problem we address is to find for fixed frequency u/27r, the monochro- 
matic solutions of the wave equation and of Maxwell's equation which achieve 
the highest field values at a point or more generally the greatest electrical 
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energy in ball of fixed small radius. They are constrained by the energy at 
\x\ ~ oo. This leads to several variational problems; 

• Maximize the field strength at a point. 

• For fixed R, maximize the energy in a ball of radius R. 

A third problem is, 

• Find Roptimai so that the energy density is largest. 

We show that the last is degenerate, the maximum occuring at R — 0. 

There are experiments in course whose strategy is to focus a number of 
coherent high power laser beams on a small volume to achieve very high 
energy densities. The problem was proposed to me by G. Mourou because 
of his leadership role in the European Extreme Light Initiative. If by bet- 
ter focusing one can reduce the size of the incoming lasers there would be 
significant benefits. Identification of the extrema guides the deployment of 
the lasers. With the experimental design as motivation the maximization 
for the Maxwell equations is interpreted as maximization of the energy in 
the electric field E, ignoring the magnetic contribution. Including the mag- 
netic contribution creates an analogous problem amenable to the techniques 
introduced here. 

As natural as these questions appear, I have been unable to find previous 
work on them. 

Study solutions of the scalar wave equation and of Maxwell equations, 

v tt - Av = 0, Et = curlS, B t = -cur\E, div E = div B = , 

for spatial dimensions d > 2. Units are chosen so that the propagation speed 
is 1. 

Definition 1.1 A solution v of the wave equation is monochromatic if it is 
of the form 

v = if>(t)u(x), with i/j" + u?i> = 0, u > 0. (1.1) 
Monochromatic solutions of Maxwell's equation are those of the form 

ip(t) (E(x) , B(x)) , with f + = 0, to > 0. (1.2) 
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They are generated by 

e ±iut u{x), and e ±wt (E(x) , B(xj) . (1.3) 

Scaling t, x — > ujt,ujx reduces the study to the case u — 1. In that case, the 
reduced wave equations are satisfied, 

(A + l)u(x) = 0, (A + l)E(x) = (A + 1)50) = 0. (1.4) 

Notation. The absolute value sign \ \ is used to denote the modulus of 
complex numbers, the length of vectors in C d , surface area, and, volume. 
Examples: l^" 1 ] and \B R (0)\. 

Example 1.2 The plane waves e*' ±w * + £ a: ) with |£| = tu is a monochromatic 
solution of the wave equation. Its period and wavelength are equal to 2tt/u. 
For Maxwell's equations the analogue is E = e l ( ±ult +i x ) e with e E C d satisfy- 
ing £ • e = to guarantee the divergence free condition. 

The solutions that interest us tend to zero as \x\ — > oo. 

Example 1.3 When d = 3, u(x) := sin|x|/|x| is a solution of the reduced 
wave equation (see also Example \3.8\) . The corresponding solutions of the 
wave equation is 



sin 


x\ 




X 





, e i(±t+\x\) e i(±t-\x\) 

2 V |x| \x\ 

For the plus sign, the first term represents an incoming spherical wave and 
the second outgoing. To create such a solution it suffices to generate the 
incoming wave. The outgoing wave with the change of sign is then generated 
by that wave after it focuses at the origin. 

Example 1.4 Finite energy solutions of Maxwell's equations are those for 
which J Rd \E\ 2 + \B\ 2 dx < oo. They satisfy 

\/R>0, lim / \E(t,x)\ 2 + \B(t,x)\ 2 dx = 0. 

J\x\<R 

Therefore, the solution (E, B) = is the only monochromatic solution of 
finite energy. 



3 



The solutions (E(x), B(x)) that tend to zero as x — > oo define tempered 
distributions on M. d . When (E(x), B(x)) is a tempered solution of the reduced 
wave equation, the Fourier Transforms satisfy 

(i-KI 2 )%) = (i-KI 2 )%) = 0. 

Therefore the support of E is contained in the unit sphere S^ 1 := {|£| = l}. 
Since 1 — |£| 2 has nonvanishing gradient on this set it follows that the value 
of E on a test function is determined by the restriction of if} to S 1 . 
Therefore there is a distribution e G X>'({|£| = 1}) so that 

E(x) := / e^e(Oda, (1.5) 

yiei=i 

where we use the usual abuse of notation indicating as an integral the pairing 
of the distribution e with the test function e^L,^ Conversely, every such 
expression is a tempered vector valued solution of the reduced wave equation. 

If e is smooth, the principal of stationary phase (see §2.21) shows that as 
\x\ — > oo, 

E(x) = ^0j- 2 (e"*N e((-x/|x|) + e^" 1 ^ e <W e(x/\x\) + 0(l/\x\) 
The field in 0(|r|~^ -1 ^ 2 ). In particular, 



;i.6) 



sup R' 1 / \E(x)\ 2 da < oo , and. (1.7) 

#>1 J\x\<R 

lim / \E(x)\ 2 dx = c d |e(£)| 2 da . (1.8) 

R -*°° JR<\x\<2R •'151=1 

For a field defined by a distribution e, (11.7p holds if and only if e e L 2 (S" i_1 ). 
In that case (II. 8p holds and stationary phase approximation holds in an L 2 
sense. This is the class of solutions of the reduced wave equation that we 
study. Equation fjl.Sp shows that HeH^^d-i) is a natural measure of the 
strength of the field at infinity. 

The divergence free condition in Maxwell's equations is satisfied if and only 
if £ ■ e(£) = on S^ -1 . In that case, the solutions e ±lt E(x) of the time 
dependent equation are linear combinations of the plane waves in Example 
PI Denote by H the closed subspace of e e L 2 (5 d " 1 ; C d ) with f • e = 0. 
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For e G H and x = r£ with |£| = 1 and r >> 1, the solution e lt E(x) of 
Maxwell's equations satisfies 

rW*E(x) « -^fe^e(-e) + e^" 1 ^ 4 e(£)V 

In practice, the incoming wave 

is generated at large r and the monochromatic solution is observed for t » 1. 
The phase factor e™'' 1 " 1 '/ 4 correponds to the phase shift from the focusing 
at the origin. 

The first two variational problems for Maxwell's equations seek to maximize 
Ji(e) := \E(0)\ 2 , and J 2 (e) := / \E(x)\ 2 dx . 

J\x\<R 

among e G H with J^-i |e(£)| 2 <icr = 1- 

Theorems 13.11 and 13.21 compute the maxima of J\ in the scalar and electro- 
magnetic cases. The maximum in the scalar case and also the vector case 
wiithout divergence free condition is IS* 01-1 !. It is attained when and only 
when e is constant. For the electromagnetic case, £ -e(£) = so the constant 
densities are excluded. The maximum is achieved at multiplies and rotations 
of the field £(£) from the next definition. 

Definition 1.5 For £ G S^ 1 denote by the projection of the vector 
(1,0,..., 0) orthogonal to £ ; 

m := (i,o,. ..,o) - (£-(i,o,...,o))e = (i,o,. ..,o) - ae- (1.9) 

i is a vector field whose integral curves are the lines of longitude connecting 
the pole (— 1, 0, • • • ,0) to opposite pole (1, 0, • • • ,0). The maximum value of 
J\ for electromagnetic waves is smaller by (d — l)/d than the extremum in the 
scalar case. The same functions also solve the J2 problem when R is not too 
large. The study of J\ is reduced to an application of the Cauchy-Schwartz 
inequality. 

In £jH the maximization of J2 is transformed to a problem in spectral theory. 
Maximizing J2 is equivalent to finding the norm of an operator. In the 
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scalar case we call the operator L. Finding the norm is equivalent to finding 
the spectral radius of the self adjoint operator L*L. The operator L*L is 
compact and rotation invariant on L 2 (S' a! ~ 1 ). Its spectral theory is reduced 
by the spaces of spherical harmonics of order k. On the space of spherical 
harmonics of degree k, L*L is multiplication by a constant A^k(R) computed 
exactly in terms of Bessel functions in Theorem 15.21 Theorem 7.1 shows that 
in the scalar case A^ j0 (-R) is the largest for R < it/2. 

For the Maxwell problem, the corresponding operator is denoted L* M L M . 
We do not know all its eigenvalues. However two explicit eigenvalues are 
(2/3)(A d) o(-R) - Ad,2(-R)), and, A di i(R). When R < vr/2, Theorem 7.3 proves 
that the they are the largest and second largest eigenvalues. The proof uses 
the minimax principal. The eigenfunctions for the largest eigenvalue are 
rotates of multiplies of I. 

For R > it/2 we derive in §8 rigorous sufficient conditions guaranteeing that 
the same functions provide the extremizers. The conditions involve the Ad t k- 
To verify them we evaluate the integrals defining the A^,k approximately. 
By such evaluations we show that for d = 3 and R < 2.5 the solutions 
maximizing J\ also maximize J 2- 

The energy density is equal to the largest eigenvalue divided by |-Br(0) |. In 
the range d = 3, R < 2.5 in both the scalar and electromagnetic case this 
quantity is a decreasing functions of R. This shows that the third problem 
at the start, of finding the radius with highest energy density is solved by 
R = 0. However the graph is fairly flat. The density dips to about 1/2 its 
maximum at about R = 2 which is about a third the wavelength. 

For focusing of electromagnetic waves to a ball of radius R no larger than one 
third of a wavelength the optimal strategy is to choose e(£) a multiple of a 
rotate of The extremizing electric fields are as polarized as a divergence 
free field can be. When d = 3 formula f ll.6p and Example 13.101 show that the 
far field for this choice is equal up to rotations by 



This field is cylindrically symmetric with axis of symmetry along the xi-axis. 
The restriction of £ to the unit sphere is cylindrically symmetric given by 
rotating the following figure about the horizontal axis. 

For the problem of focusing a family of lasers, this suggests using linearly 
polarized sources concentrated near vertical equator and sparse near the poles 
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on the horizontal axis. In contrast, for scalar waves one should distribute 
sources as uniformly as possible. 
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and electromagnetic cases respectively. J. Szeftel, G. Allaire, C. Sogge, P. 
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Nord, and, Universita di Pisa. Sogge and Gerard were guests of the Centro 
De Giorgi. Schotland and I were both guests of the IMA and Michigan State 
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2 Monochromatic waves. 

2.1 Electromagnetic waves and their transforms. 
Proposition 2.1 i. E given by (11.51) satisfies div E = if and only if 



ii. If a monochromatic solution of the Maxwell equations has electric field 
given by (11.31) with u — 1 and E is given by (II. 5p then the magnetic field is 
equal to e %i B[x) with 



e(0 " t 







on 



{iei = i}. 



(2.1) 




(2.2) 



Proof. Differentiating ( II. 5p yields, 

divE = [ e'^ i£-e(£)da, curl£ = / e ix ^ i£ A e(f) da . 
The first formula proves i. 

The Maxwell equations together with (11.31) yield 

-curlE = B t = %B. 
Therefore, the second formula proves ii. I 

Remark 2.2 The condition (12. ip asserts that e(£) is tangent to the unit 
sphere. Brouwer's Theorem asserts that if £ i-> e(£) zs continuous then there 
must be a £ where e(£) = 0. 

Example 2.3 . If d = 3 and E is given (II. 5p e(£) = i/ien on 
1*1 = 1, 

^A£(0 = £ A ((1,0,0) - 6^) = £ A (1,0,0) = (0,6,-6) 

is the tangent field to latitude lines winding around the x±-axis. Since this is 
an odd function, the magnetic field vanishes at the origin, B(0) = 0. 

2.2 The competing solutions. 

First verify the stationary phase formula ( 11. 6ft from the appendix. Consider 
E given by ( 11. 5ft with e G C°°({|£| = 1}). For x large, the integral ( II. 5p has 
two stationary points, £ = ±x/|x|. At £ = x/\x\ parameterize the surface by 
coordinates in the tangent plane at x to find that the phase x£ has a strict 
maximum equal to \x\ and hessian equal to the — /(d-i)x(d-i)- At £ = —x/\x\ 
the phase has a minimum with value — \x\ and hessian equal to the identity. 
The stationary phase method yields (II. 6p . The energy in the electric field 
satisfies ([LTD and Oil . 

Theorem 2.4 Suppose that E G 5'(IR d ) is a tempered solution of the reduced 
wave equation given by ( 11.51) with e G T>\S d ~ l ). Then, ( j 1.7ft ZioZds if and only 
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ife G L 2 (S d 1 ; C d ). In that case fll.Sp holds. In addition the stationary phase 
approximation holds in the sense that as R —> oo, 



I 



R<\x\<2R 



E (^)- 1(^/2 ( e " N e((-x/|x|)+e^- 1 )/ 4 e^l e(x/|x|)) f dx = o(R) 



Proof. The first two assertions are consequences of Hormander [lj Theorems 
7.1.27 and 7.1.28. The Theorem 7.1.28 also implies that 

1 f . |2 f |2 

sup — / dx < c{d) I | e (0| da. 

R>1 R JR<\x\<2R [€|=1 

This estimate shows that to prove the third assertion it suffices to prove it 
for the dense set of e 6 C°°(5' <i_1 ). In that case the result is a consequence 
of the stationary phase formula (jl.6p . I 

Definition 2.5 H is the closed subspace of e G L 2 (S' d_1 ; C d ) consisting of 
e(£) so that £ ■ e(£) = 0. Denote by U the orthogonal projection of e G 
L 2 (S d ~ l ; C d ) on H. 

The next example explains a connection between the solutions of the reduced 
equation that we consider and those satisfying the Sommerfeld radiation 
conditions. 

Example 2.6 If g G £'(1^ : C d ) is a distribution with compact support, then 
there are unique solutions of the reduced wave equation 

(A + l)E out = g, (resp. (A + l)E in = g) 

satsifying the outgoing (resp. incoming) radiation conditions. The difference 
E := E out — Ei n is a solution of the homogeneous reduced wave equation. The 
field F := e lt E(x) is the unique solution of the initial value problem 

OF = 0, F| t=Q = g, F t \ f=0 = ig. 

When g G L 2 the formula E(x)5{r — 1) = c e~ trt Fdt together with the 
solution formula for the Cauchy problem imply that (11.71) holds (or see 
Theorem 14-3.4 showing that both incoming and outgoing fiels satisfy ( 11. 71) ). 
More generally, the Fourier transforms in time of solutions of Maxwell's 
equations with compactly supported divergence free square integrable initial 
data yield examples of monochromatic solutions in our class. 
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2.3 Spherical symmetry is impossible. 



It is natural to think that focusing is maximized if waves come in equally in 
all directions. For the scalar wave equation that is the case. However, such 
waves do not exist for Maxwell's equations. Whatever is the definition of 
spherical symmetry, such a field must satisfy the hypotheses of the following 
theorem. 

Theorem 2.7 // E(x) G C 1 ^) satisifies div E = and for x ^ the 
angular part 

V fI/ \ x \ 

has length that depends only on \x\ then E is identically equal to zero. 

Proof. The restriction of the angular part of E to each sphere \x\ = r is 
a C 1 vector field tangent to the sphere and of constant length. Brouwer's 
Theorem asserts that there is a point x on the sphere where the tangent 
vector field vanishes. Therefore the constant length is equal to zero and E is 
radial. 

Therefore in x ^ 0, 

E(x) = 4>(\x\) x . 
Since E G C 1 it follows that G C\{\x\ > 0}). 
Compute for those x, 

div E = 4>divx + (V x (p)-x = d(p + r (fi r . 

Therefore in x ^ 0, r <ft r + d<fi = so <ft — cr~ d . 

Since E is continuous at the origin it follows that c = so E = in x ^ 0. 
By continuity, E vanishes identically. I 



3 Maximum field strengths. 

We solve the variational problems associated to the functional J\ to yield 
sharp pointwise bounds on monochromatic waves. The fact that the bounds 
for electromagnetic fields are smaller shows that focusing effects are weaker. 
The extremizing fields are first characterized by their Fourier Transforms. 
Explicit formulas in x-space are given in §3.41 
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3.1 Scalar waves. 
Theorem 3.1 If 

[ <**f(t)da, f G L 2 (S 
J\e\=i 



111 x } 



and x G R , then 



lulx) 



lL 2 (5 d -!) 

with equality achieved if an only if f is a scalar multiple of e~ 1 -^ . 

Proof. The quantity to maximize is the L 2 (S d ~ 1 ) scalar product of / with 
e -iz£ rj^g resu it j s exactly the Cauchy-Schwartz inequality. I 

3.2 Electromagnetic waves. 

From Definition II .5[ is tangent to the longitude lines on the unit sphere 
connecting the pole (—1, 0, . . . , 0) to the pole (1, 0, . . . , 0). It is the gradient 
of the restriction of the function £ x to the unit sphere. 

Theorem 3.2 If d > 2 and 

E(x) = [ e ixi e(0 da, e G H . 



Then 



'd-1 \ 1/2 



Equality holds if and only if e is equal to a constant mulitple of a rotate of 

no- 

Proof. By homogeneity it suffices to consider ||e||£2(gd-i) = 1. Rotation 
and multiplication by a complex number of modulus one reduces to the case 
£7(0) = |£(0)|(1, 0, . . . , 0) and |£7(0)| = / ei(£) da. Need to study 



sup 



ei (0da : e-e(0 = 0, J _MQ\\ 2 da = l} 
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The quantity to be maximized is 

/ ei(£)da = (e, (1,0. ■ ■ • ,0)) L2fsd _ M 
J\C\=i 

The constant function, (1, 0. ■ ■ • ,0) does not belong to the subspace H. The 
projection theorem shows that the quantity is maximized for e proportional 
to the projection of (1, 0. ■ • • , 0) on H. Equivalently, using (jl.9p together 
with e • £ = yields 

ei = e-(l,0,...,0) = e.(*(£) + CxC) = e-£. 



L 



ei(0do- = / e(0-m da (3.2) 



|5|=i ■/ [€l=i 

which is equal to the H scalar product of e and £. Since one has the orthog- 
onal decomposition 

(1,0,..., 0) = *(0 + £i£, one has, 1 = |£(£)| 2 + & • 
The Cauchy-Schwartz inequality shows that the quantity (13. 2p is 

< ||e|| £ a(5d-i) pHia^d-i) = WeWrfi^i^J (l-^do^J . (3.3) 

The extremum is attained uniquely when e = z£/\\£\\ with \z\ = 1. 
To evaluate the integral on the right of (13.31) compute, 



da 



\s- 



d-l\ 



|5|=1 J |f | = 1 UJ\(\; 

Therefore 



d 



;i-e?)do- = is"- 1 ! - ^— ^ = \s d - x \ 



'ICI=l 

Together with (I3.3P this proves (13. ip . 



Remark 3.3 // one constrains e to /jave support in a subset Q then with x 
denoting the characteristic function of Q, 

f e(0 • *(£) do- = / e(0 • £(0x(0 da 

and Ei is maximized by the choice e = £(£)x(0- ^ n ^ e extreme light initia- 
tive Q is a small number of disks distributed around the equator xi — 0. 
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3.3 Derivative bounds. 

Corollary 3.4 If d > 2 and E satisfies 

E(x) = [ e Lxi e(f) da, e e H, 

then for all a E N d and x £ M d , 

( d — 1 \ 1//2 

Proof. The case a = follows from Theorem 13.21 applied to 

E(x) := E{x + x) = [ e** e** e(£) da := / e^e(0^- 

7iei=i 7|ei=i 

Compute for |a| > 

d a x E = d% [ e**e(£)d(T= [ e ix Z {i£) a e{£) da , 

J\i\=i J\e\=i 

which is of the same form as E with density (i£) a e(£) orthogonal to £. Since 
|£j| < 1 it follows that < 1 so || (iO^IU 2 ^- 1 ) — ll e llL 2 (5 d - 1 )- Therefore 
the general case follows from the case a = 0. I 

Remark 3.5 Derivative bounds for the scalar case are derived in the same 
way. They lack the factor (d — l)/d. 

3.4 Formulas for the extremizing fields. 

The electric field corresponding to the extremizing density I is explicitly 
calculated. The computation relies on relations between Bessel functions, 
spherical harmonics, and, the Fourier Transform. These relations are needed 
to analyse J 2 . 

Start from identities in Stein- Weiss [2]. Their Fourier transform is defined 
on page 2, 

J f(x) e~ t2nx ^ dx, n.b. the 2n in the exponent. 
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We will not follow this convention, so adapt their identities. The Bessel 
function of order k is (page 153), 



r[(2fc + i)/2]r(i/2) 

Theorem 3.10 (page 158) is the following. 



j 1 e lts (1 - s 2 )^- 1 )/ 2 ds, -l/2<keR. 

(3.5) 



Theorem 3.6 If x E W d , f = f(\x\)P(x) E L x (]R d ) with P a homogeneous 
harmonic polynomial of degree k, then j f(x)e~ 2mx ^dx = _F(|£|)P(£) with 

POO 

F(r) = 2m~ k r^ d+2k - 2 ^ 2 / f(s) J {d+2k _ 2)/2 (2nrs) s {d+2k)/2 ds . 



This theorem is equivalent, by scaling and linear combination, to the same 
formula with / = 5(r — 1). That case is the identity, 

/ e~ i2 ^ P(x) da = 2m- k \^ d+2k ~ 2)/2 J {d+2k - 2)/2 {2-n\i\) P(0- (3-6) 
.71*1=1 

Remark 3.7 i. For |f | ->■ oo ; J(|£|) = O^' 1 ' 2 ), and P{£) = 0(\£\ k \) so 
the right hand side is 0(|^| - ^ 2 ^ 2-1 / 2 ) = 0(|^| - ^ -1 ^ 2 ) as required by the 
principle of stationary phase. 

ii. For |f| -»■ 0, J((d- 2 fc-2)/2(|C|) = 0(|^|( d - 2fc - 2 )/ 2 ) so the right hand side 
of (I3.6P is 0(|£| fc ). The higher the order of P the smaller is the Fourier 
transform near the origin. 

To adapt to the Fourier transform without the 27r in the exponent, use the 
substitution 77 = 27r£, \r]\ = 2n\t;\ to find, 

e~^P{x)da = 27rr fc (|r ? |/27r)-( rf+2fc - 2 )/ 2 J (d+2fc _ 2)/2 (|r ? |)P(r / /27r). 

31=1 

Using the homogeniety of P yields 

= {2^f- k i- k {\r l \/2^ d+2k - 2 ^ 2 J (d+2fe _ 2)/2 (M) P(v). 
The exponent of 2ir is equal to d/2 yielding, 

e** P(x) da = (2rr) d / 2 r k \ v \-( d + 2k ^/ 2 J {d+2k - 2)/2 (\ V \) P(r,). (3.7) 



X =1 
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Since, \r]\ k P(r]) = P(r)/\r]\) (13.71) equivalent to, 

/ e-^ P{x) da={2^) d lH- k \r 1 \-^l 2 J {d+2k _ 2 y 2 {\r 1 \) P( V /\v\)- (3.8) 
71*1=1 

The change of variable 77 1— > —77 yields, 




P(x) da = (2vr) d / 2 (-z)~* \ V \-^ 2 J (d+2k - 2)/2 (\ V \) P(v/\v\)- 



(3.9) 

Finally interchange the role of x and rj to find, 

/ e^ P( V ) da = (2n) d / 2 (-i)~ k \x\-^/ 2 J {d+2k - 2)/2 {\x\) P{x/\x\). 
J\n\=x 

(3.10) 

Example 3.8 The second most interesting example is the extremizing field 
for the scalar case when d = 3. In that case P = constant and there is a 
short derivation. The function u(x) := j^ =l e lx ^da is a radial solution of 
(A + l)u = 0. In x 7^ these are spanned for d = 3 by e ±ir jr. Smoothness at 
the origin forces u = Asinr/r. Sinceu(0) = \S d ~ 1 \ it follows that A = \S |. 

The most interesting case for us is d = 3 and the extremizing field E with 
e(£) = £(£). Since £ is not a spherical harmonic, the preceding result does 
not apply directly. To find the exact electric field, decompose i in spherical 
harmonics. 

Lemma 3.9 The spherical harmonic expansion of the restriction of £((,) to 
the unit sphere S C M d is 




on |e| = l. (3.11) 



Proof of lemma. When |£| = 1, 

£(0 = (1, 0, . . . , 0) - fif = (1, ■ ■ ■ , -ted) - (& 0, . . . , 0) . (3.12) 

The first summand has coordinates that are spherical harmonics. 
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Decompose 

to find the expansion in spherical harmonics of the restriction of ^ to the 
unit sphere £f + • — h £j = 1, 

£i = 3+E^ on e? + - + e2 = l- (3-13) 

i=2 

Using fl37T3|) in f[3~T2j) proves (13TTT]) . I 

Example 3.10 The stationary phase formula (11. 6j) applied to the extremiz- 
ing e = £(£) which is an even function yields for d = 3 



a/2tt V|x|/ \x\ 

This is the sum of incoming and outgoing waves with spherical wave fronts 
and each with profile on large spheres proportional to £(x/\x\). The desired 
incoming wave is such an [-wave. 



4 Equivalent selfadjoint eigenvalue problems. 

The section introduces eigenvalue problems equivalent to the maximization 
of J 2 . 



4.1 The eigenvalue problem for focusing scalar waves. 

Definition 4.1 For R > define the compact linear operator L : L 2 {S d ^ 1 ) — > 
L 2 (B R (0)) by 

(Lf)(x) := /' e«/(0<fa. 

J\i\=i 

The operator L commutes with rotations. The adjoint L* maps L 2 (Br(0)) — > 
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Proposition 4.2 The following four problems are equivalent. 

i. Maximize the functional J 2 on scalar monochromatic waves. 

ii. Find f G L 2 (S d ^ 1 ) with ||/||.L2(,s d - 1 ) = I so that ||-^/||b h (o) is largest. 

iii. Find the norm of L. 

iv. Find the largest eigenvalue of the positive compact self adjoint operator 
L*L on L 2 (S d - r ). 

Proof. The equivalence of the first three follows from the definitions. The 
equivalence with the third follows from the identity 

W\\h { B R ) = (Lf,Lf) LHBR) = (L*L/,/) L2(5d - 1} . I 

Definition 4.3 Define a rank one operator 

L'iS^-X) 3 f -> L f := / f(Oda E C. 

J\Z\=i 

Remark 4.4 i. The problem of maximizing J\ for scalar waves is equivalent 
to finding the norm of L and also finding the largest eigenvalue of LqL q . 

ii. The same formula defines an operator from L 2 (S d ~ v ) — > L 2 (Br(0)) map- 
ping f to a constant function. With only small risk of confusion we use the 
same symbol L Q for that operator too. 

Definition 4.5 The vector valued version of L and L are defined by 
(Le)(x) := [ e ix ^e^)da, e e L 2 (S d ~ 1 ; C d ) , 

L e := / e(Z)da, e e L 2 (S d ~ 1 ; C d ) . 
J\Z\=i 

Denote by L M and L M the restriction of~L and L to H. 

Remark 4.6 i. The problem of maximizing the functional J\ for monochro- 
matic electromagnetic waves is equivalent to finding the largest eigenvalue of 

ii. The problem of maximizing the functional J2 for monochromatic electro- 
magnetic waves is equivalent to finding the largest eigenvalue o/L^Lm. 

iii. The operator LTI is equal to Lm on H and equal to zero on H ± . Therefore 
(Ln)*(LIl) is equal to L* M L M on H and equal to zero on H^. 
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5 Exact eigenvalue computations. 



5.1 The operators L^Lq, LqLq and Lq M Lq 7 m- 

Theorem 5.1 i. The spectrum of LqL contains one nonzero eigenvalue, 
(S^ -1 ), with multiplicity one. The eigenvectors are the constant functions. 

ii. The spectrum of the operator LqL contains one nonzero eigenvalue, 
IS^ 1 ], with multiplicity d. The eigenvectors are the C d -valued constant func- 
tions. 

iii. The spectrum of Lq m L 0i a/ contains one nonzero eigenvalue, IS^ -1 ^ — 
l)/d with multiplicity d. The corresponding eigenspace consists of £(£,) and 
functions obtained by rotation and scalar multiplication. 

Proof iii. Suppose that / is an eigenfunction in H so that the norm of J f da 
is maximal. Rotating / yields a function in H with the same ||L /|| and with 
Lq/ parallel to (1, 0, . . . , 0). Therefore maximizing Lq/ and maximizing 

L /-(1,0,... ,0) 

yield the same extreme value. 
Since / G H, 

f fida = y /- (1,0, ...,0)) da = J f ■ £ da . 

The extreme value is attained for / parallel to t. This shows that i is an 
eigenfunction corresponding to the largest eigenvalue. 

Rotating and taking scalar multiples yields a complex eigenspace of dimen- 
sion d. Since rank L = d these are all the eigenfunctions. 

If ||/||i2(sd-i) = 1, then the maximization of Jo shows that |L ,m/| 2 = 
\S d ~ 1 \(d — l)/d proving the formula for the eigenvalue. I 

5.2 The eigenfunctions and eigenvalues of L*L and L*L. 

Theorem 5.2 In dimension d, the spherical harmonics of order k are eigen- 
functions of L*L with eigenvalue 

A d , k (R) := (2 7 r) d / 2 |^ 1 | fr [J {d+2k - 2)/2 (r)} 2 dr . (5.1) 

Jo 
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Remark 5.3 From the point of view of focusing of energy into balls, all 
spherical harmonics of the same order are equivalent. 

Proof. If P is a homogeneous harmonic polynomial of degree k formula 
(I3.10p shows that 

{LP){x) = <f> d)k {\x\) P{x), 

defining the function </>. 

The operator L* is an integral operator from L 2 (Br) — > L 2 (S' d ~ 1 ) with kernel 
e~ lx ^. Therefore 

L*L(P) = [ e-^<t> d , k {\x\)P{x)dx. 

J\x\<R 

Introduce polar coordinates x = ry with \y\ = 1 to find 

jd-l\ f f d-l e -iryt, 



L*L(P) = l^- 1 ! / / r d ~ x <f> d , k (r) r* P(y)da(y)dr. 

Jo J\y\ = l 

Formula f)3.10p shows that 

e-^ P(y) da(y) = <f> d , k {r) P(-r £) = {-r) k <j> dik (r) P(0 . 



Ii/I=i 
Therefore 

L*L(P) = A-l) fe |^V d+2fe_1 ^ lfe (r) 2 drP(e). 
Jo 

This proves that the harmonic polynomial are eigenvectors with eigenvalue 
depending only on d and k. 

The formula for the eigenvalue follows on noting that the eigenvalue is equal 
to the square of the L 2 (B R (0)) norm of 

F(x):=[ e ixr > P(ri) da , [ \P(rj) | 2 da = 1. (5.2) 
J\ v \=i J\ v \=i 

Using polar coordinates and f)3.10p yields, 



\ F \\l R (o) = (2vr) d/2 |^- 1 | J^r-^ [<W- 2)/2 (r)] 2 r"" 1 dr 

= ^nf 2 ^ f R r [J {d+2k . 2)/2 {r)] 2 dr . 
Jo 
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(5.3) 



This proves (15. ip . 
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The spectral decomposition of L is nearly identical to that of L. The next 
result is elementary. 

Corollary 5.4 The eigenvalues of L*L are the same as the eigenvalues of 
L*L. The eigenspaces consists of vector valued functions each of whose com- 
ponents belongs to the corresponding eigenspace of L*L. 

5.3 Some eigenf unctions and eigenvalues of h* M h M . 

The situation for L* M L M is more subtle. Our first two results show that there 
are eigenfunctions intimately related to the eigenvalus A d l (i?) and A d (-R)- 

Theorem 5.5 The d dimensional space of functions e(£) := £ A £ with ( G 
C d \ are eigenfunctions of~L* M ~L M . The eigenvalue is A d l (i?). 

Remark 5.6 These e(£) are the vector valued spherical harmonics of degree 
1 that belong to H ; that is, that satisfy £ • e(£) = 0. 

Proof. Follows from Le = A d l (R) e and e G H. I 

Though the constant functions which are eigenvectors of L do not belong to 
H, their projection on H yield eigenvectors of L* M ~L M . 

Theorem 5.7 The d dimensional space consisting of scalar multiples of ro- 
tates of consists of eigenfunctions of~L* M \j M with eigenvalue equal to 

(A dfi (R) - A d , 2 (R)) ^=i. (5.4) 

Proof of theorem. Use (13. lip . Since the spherical harmonics are eigen- 
functions of L*L one has, suppressing the R dependence of A, 

rl 1 d — £ 2 

L*L t = (A d ,o— p - A d , 2 ~ L ^ L > 66 , • • • , -Ad >2 Ud) ■ 

i=2 
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Multiply (IXTTj) by A d>2 to find on + h £j = 1, 

A d , 2 ^ = (A d , 2 ^ - A d , 2 ^ , -A d)2 &6 , ■ ■ • , -A dj2 a^d) • 

V i=2 ' 

Subtract from the preceding identity to find, 

L*L£ = (A d , -A d , 2 )^(l,0,...,0). 

Projecting perpendicular to £ using II (1, 0, • • ■ ,0) = £ yields 

ni/LIK = ITI/L£ = (A dfi -A d ^^± i. 

This proves that I is an eigenfunction of (LIT)* (LIT) with eigenvalue (A^o — 
A d , 2 )(d-l)/d. 

Remark I4.6l iii shows that it is an eigenfunction of L^ f L M with the same 
eigenvalue. 

By rotation invariance the same is true of all scalar multiples of rotates of £. 
They form a d dimensional vector space, spanned by the projections tangent 
to the unit sphere of the unit vectors along the coordinate axes. I 

As in Theorem 15.51 if one defines to consist of spherical harmonics of 
degree k that belong to H, then are orthogonal eigenspaces of L* M ~L M 
with eigenvalue Ad,k(R)- 

Example 5.8 In M. 2 the homogeneous M 2 valued polynomials of degree two 
whose radial components vanish are spanned by (— Xix 2 ,x 2 ) and (x 2 ,, —X1X2). 
There are no harmonic functions is their span proving that when d = 2, 
H 2 = 0. It is clear that H = 0. 

Though there are a substantial number of eigenvectors of L* M ~L M accounted 
for by the H fc they are far from the whole story. 

6 Spectral asymptotics. 

6.1 Behavior of the A^(i^). 

Proposition 6.1 i. As R ->■ 0, A d>k (R) = 0(R d+2k ). 



21 



ii. As R -> 0, A d)0 (R) = {S^l \B R {0)\(1 + 0{R)). 

iii. linifc^oo A^k{R) = uniformly on compact sets of R. 

Proof, i. Formula (13.51) shows that Jh(t) = 0(t k ) as t — > 0. Assertion i then 
follows from (15. ip . 

ii. By definition, Ad t o(R) is the square of the L 2 (Br) norm of Jj^ =1 e lx ^ /(£) der 
for / a constant function of norm 1. Take / = \S d ~ l \~ l l 2 . 

For R small e** = 1 + O(R) so 

/ e ix ^f(Oda= [ (l + 0(R))\S dr - 1 \- 1/2 da=\S dr ' 1 \ 1/2 {l + 0(R)). 
J\£\=l J\£\ =1 

Squaring and integrating over Br proves ii. 

iii. As k — > oo, the pref actors in the formula for tend to zero uniformly 
since T([(2k + l)/2] dominates. 

The integral in the definition tends to zero uniformly by Lebesgue's Domi- 
nated Convergence Theorem. I 



6.2 Small R asymptotics of the largest eigenvalues. 

Each of the operators L, L, and ~L>m has integral kernel e lx ^. They differ in 
the Hilbert space on which they act. For x small, e tx ^ ~ 1 showing that for R 
small the three operators are approximated by L , L , and L 0) m respectively. 
We know the exact spectral decomposition of the approximating operators. 
Each has exactly one nonzero eigenvalue. In performing the approximation 
some care must be exercised since the operator to be approximated has norm 
0(R d/2 ) tending to zero as R — > 0. 

Proposition 6.2 i. Each of the operators L, L, and ~Lm has norm no larger 
than(\B R (0)\\S d - 1 \) 1 / 2 . 

ii Each of the differences L — Lq, L — Lo, and ~Lm — Lo,m has norm no larger 
than 

(d + 2)V2 ■ (6-1) 

Proof, i. Treat the case of L. For ||e|| = 1, the Cauchy-Schwartz inequality 
implies that for each x, ||Le(x)|| 2 :d < \S |. Integrating over the ball of 
radius R proves i. 
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ii. The Cauchy-Schwarz inequality estimates the difference by 



:L-L n )e 



I (e** - 1) e(0 da < [ \x\ |e| da 



For e of norm one this yields 



- L )e||^2( Bfl ( )) < \S d x \ 



\x\ dx 



B T 



R 



\S d ~^\ I I r 2 r* -1 dr da(u) 

'\uj\=1 JO 



Ud+2 
d-l 1 2 JX 



d+2 



completing the proof. I 

Theorem 6.3 i. For each d there is an Ri{d) > so thatforO < R < -Rl(^) 
the eigenvalue A^o(-R) is the largest eigenvalue of L*L. It has multiplicity 
one. The eigenfunctions are constants. 

ii. For < R < R^d) the eigenvalue A di0 (R) is the largest eigenvalue of 
L*L. It has multiplicity d. The eigenfunctions are constant vectors. 

iii. For each d there is an -Rm(^) > so that for < R < Ru(d) the 
eigenvalue 



(A d ,o(R) - A d ,2(R)) 



d- 1 
d 



(6.2) 



is the largest eigenvalue o/I/^L M . It has multiplicity d. The eigenfunctions 
are rotates of constant multiples of I. 

iv. In all three cases, the other eigenvalues are 0(R d+1 ). 



Proof. We prove iii and iv for the operator L* M ~L M . Proposition 16.21 im- 
plies that the compact self adjoint operators L* M L M and Lq a/ L m differ by 
0(R d+1 ) in norm. 

Part iii of Theorem 15.11 shows that the spectrum of Lg M L M contains one 
positive eigenvalue, A + := |-Br(0)| |S' a! ~ 1 |((i— l)/d. The factor |_B^(0)| arises 
because L 0iM in the present context is viewed as an operator with values in 
the functions on -Br(O) C M d . The eigenfunctions are scalar multiplies of 
rotates of I. The rest of the spectrum is the eigenvalue 0. 
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It follows that the spectrum of ~L* M ~L M lies in the union of disks of radius 
0(R d+1 ) centered at zero and A + . For R small these disks are disjoint and 
the eigenspace associated to the disk about A + has dimension d. 

Theorem 15.71 shows that the eigenf unctions of Lq m L m with eigenvalue A+ 
are eigenfunctions of L A/ L M . The eigenvalue is given by (15. 4ft . 

It follows that for R small, the scalar multiples of rotates of £ is an eigenspace 
of 

La/La/ of dimension d and eigenvalue in the disk about A+. This completes 
the proof of iii. 

The fact that the other eigenvalues lie in a disk of radius 0(R d+1 ) centered 
at the origin proves iv . 

The proofs for the operators L and L are similar. I 

7 Largest eigenvalues for R < n/2. 

7.1 Monotonicity of Ad,k(R) in k. 

Recall that the wavelength is equal to 2n. 

Theorem 7.1 For < R < tt/2, Ad t k(R) is strictly monotonically decreas- 
ing in k = 0, 1, 2, ... . In particular the largest eigenvalue of L*L and L*L is 
A-d,o{R)- The corresponding eigenfunctions are constant scalar and constant 
vector functions respectively. 

Proof. Write, 

- mrmwm /' cos(rs) (1 - s2)< " 1,/2 da - (7A) 

For < r < tt/2 the cosine factor in the integral is positive. Since (1 — 
s 2^(2fc-i)/2 j g cl ecreas i n g j n k f or s £ [o ; 1] ; the integral is decreasing in k. 

F[(2k + l)/2] is increasing in k. Since r < 2, (r/2) fc decreases with k. The 
proof is complete. I 

Remark 7.2 The first figure in §EU\ shows that the graphs of the functions 
A-3,o(R) o,nd and K^^R) cross close to R = it . The proof shows that this 
cannot happen for R < tt/2. 



24 



7.2 Largest eigenvalues of L^L M when d = 3, R < tt/2. 

Theorem 7.3 When d = 3 and R < tc/2 the strictly largest eigenvalue of 
L* m Lm is (A 3j0 (-R) — A 3i2 (_R))2/3. The eig en functions are the scalar multiples 
of rotates of i, and, A 31 (i?) is the next largest eigenvalue. 

The proof uses the following criterion valid for all d, R. For ease of reading, 
the R dependence of A dk (R) is often suppressed. 

Theorem 7.4 i. The eigenvalue (A di0 — A di2 )(d — l)/d of~L* M ~L M is strictly 
larger than all others only if 

(A dfi -A dt2 )(d-l)/d > Ad,!. (7.2) 

ii. If in addition to (17.21) . the two largest evalues of L*L are A^o and A d ,x, 
then the eigenvalue [A d fi — A d ^){d— l)/d of~L* M ~L M is strictly larger than the 
others. The eigenfunctions are the scalar multiples of rotates of £, and, A d)1 
is the next largest eigenvalue of~L* M ~L M . 

Remark 7.5 Equation ( 17.2 j) implies that A^o > A^i- The additional condi- 
tion in ii is that for all k > 1, A d i > A d ^. 

Proof of Theorem 17.41 i. Since A^i is also an eigenvalue of L^L M , 
necessity is clear. 

ii. Under these hypotheses Theorem 15.21 shows that A^o is the largest eigen- 
value of L*L with one dimensional eigenspace consisting of constant func- 
tions. Corollary 15.41 shows that L*L has the same largest eigenvalue with d 
dimensional eigenspace consisting of C d valued constant functions. The next 
largest eigenvalue of L*L is A^i. In particular, L*L has exactly d eigenvalues 
counting multiplicity that are greater than A^.i. 

Since Lm is the restriction of L to a closed subspace, the minmax principal 
implies that L* M L M has at most d eigenvalues counting multiplicity that are 
greater than A^i- 

Theorem 15.71 provides a d dimensional eigenspace with eigenvalue given by 
the left hand side of (17.21) and therefore greater than A^i. In particular there 
are exactly d eigenvalues counting multiplicity that are greater than A^i- 

Theorem 15.51 shows that A rfl is an eigenfunction of L^ f L A/ so it must be the 
next largest. I 
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Example 7.6 Parts i and ii of Proposition W7T\ show that the sufficient con- 
dition is satisfied for small R. This gives a second proof that for small R, £ 
is an extreme eigenfunction for L* m ~Lm- The first proof is part iii of Theorem 

Proof of Theorem 17.31 Verify the sufficient condition of Theorem I7.4l ii. 

Since R < tt/2, A 3 ^. are strictly decreasing in k. Therefore A 3 and A 31 are 
the two largest eigenvalues of L*L. 

It remains to verify (17. 2p . Use formulas (15. ip and (17. ip . Formula (15. ip with 

d = 3 and k — 0,1,2 involves Ji/ 2 , J3/2, ^5/2- 

Since the integral in (17. ip is decreasing in k it follows that 

J k+ i{r) < r r((2A; + l)/2) 



J k {r) ~ 2 r((2fc + l)/2 +1) ' 
The functional equation F(n + 1) = (n + l)r(n) yields 
J fc+ i(r) r 1 r 



J fe (r) 2 (2fc + 3)/2 2k + 3' 
Therefore, 

^3/2 (r) r J 5/2 (r) r 

— — - < -, and, < -. 

•AM?") 4 ^3/2(r) 6 

Injecting these estimates in (15. ip yields 

^iW<^, md , ^<^. (7.3) 

Therefore, 

R 2 R 2 R 2 R 2 

A 3 ,i < ^2-A 3 , , and, A 3j2 < — A 3 ,i < A 3 , , so, 

(A 3 ,o - A 3l2 ) I ~ A 3 ,i > A 3jO - ^42^2 - §") : = A 3,oW- 

The polynomial is equal to 2/3 when R = and decreases as i? in- 

creases. To verify (17.21) it suffices to show that h(Ti/2) > 0. Since 2 > ir/2, 
h(n/2) > h{2) = 43/108 > 0. I 
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8 Numerical simulations to determine largest 
eigenvalues. 



Recall that the wavelength is equal to 2ir. In this section the dimension 
d = 3. 

Theorem 15.21 Corollary I5.4[ and Theorem 17.31 allow one in favorable cases 
to find the largest eignevalues of L*L, L*L, and L* M ~L M , by evaluating the 
intergrals defining A 3t k(R) for k = 0,1,2,.... These quantities decrease 
rapidly with k so to compute the largest ones requires little work. 

8.1 Simulations for scalar waves. 

For scalar waves the eigenvalues are exactly the A 3t k(R). For R small, they 
are monotone in k so the optimal focusing is for k = 0. Our first simulation 
(performed with the aid of Matlab) computes approximately the integrals 
defining A 3j fc(i?) for R < 2n and k = 0, 1,2,3. The resulting graphs are in 
the figure on the left. The horizontal axis is R and on the vertical axis is 
plotted the integral on the right hand side of (15. II) . that is, 

A ^(R) .. A3, fc (fl) fr_ 012 o 
(27r) 3 / 2 |5 2 | 27/2^5/2' K U,1,Z,0. 

The four curves correspond to the four values of k. The graph with the 
leftmost hump is A3 j o(i2). The graph with the hump second from the left is 
A 3j i(i2) and so on. The conclusion is that A 3)0 (i?) crosses transversaly the 
graph A 3>1 (i?) just to the right of R = 3. At that point, A 31 (i?) becomes 
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the largest. On the right is a zoom showing that the crossing is suspiciously 
close to R — ir. 

The graphs of the A 3y k{R) are a little misleading since it is not the total 
energy but the energy density that is of interest. The next figure plots as a 
function of R A 3t0 / (2 7 / 2 7r 5 / 2 ^^(0) |). The small gap near R = is because 
the division by |-Br(0)| is a sensitive operation and leads to numerical errors 
in that range. 

The energy density is greatest for balls with radius close to R — 0. The 
density drops to half its maximum value at about R = 2 which is about 1/3 
of the wavelength. 

8.2 Simulations for electromagnetic waves. 

Using Theorem 17.31 one can investigate the analogous questions for Maxwell's 
equations by manipulations of the A 3k (R). 

The simulations of the preceding subsection show that for R < 3 one has 
A 3i0 (-R) > A 3) i(i2). To show that the eigenvalue corresponding to is the 
optimum it suffices to verify (J772J) . To do so one needs to verify the positivity 
of 

7/2 5/2 /2 A 3 , (R) _ 2 A3, 2 (fl) _ A 3A (R) \ 
\3\B R (0)\ 3\B R (0)\ \B R (0)\J- 

This is a linear combination of quantities computed in the preceding sub- 
section. Its graph is plotted on the left. The graph crosses from positive to 
negative near R = 2.5. The criterion is satisfied for all R to the left of this 
crossing. 
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For R < 2.5 the energy density for the optimizing monochromatic electro- 
magnetic fields associated with is equal to 

2 A 3 ,o(i?) _ 2 A 3|2 (fl) 

3 \B R (0)\ ~ 3 \B R {0)\ 

Because of the factor 2/3 it is smaller than the density in the scalar case by 
that factor. The subtraction in the formula shows that the density drops off 
more rapidly in the electromagnetic case than in the scalar case. The graph 
of 2 -7 / 2 7r~ 5,/2 times this quantity is plotted next. 

0.055 

0.06 
0.04B 

0.04 
0.035 

0.03 
0.025 

0.02 




As in the scalar case the maximal energy density occurs on balls with radius 
near zero. The intensity drops to one half of this value to the left of R = 1.9. 
This is close to the corresponding value for scalar waves, about one third of 
a wavelength. 
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